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ABSTRACT 

SN1996cr is one of the five closest SNe to explode in the past 30 years. Due to its 
fortuitous location in the Circinus Galaxy at ^ 3.7 Mpc, there is a wealth of recently 
acquired and serendipitous archival data available to piece together its evolution over 
the past decade, including a recent 485 ks Chandra HETG spectrum. In order to inter- 
pret this data, we have explored hydrodynamic simulations, followed by computations 
of simulated spectra and light curves under non-equilibrium ionization conditions, 
and directly compared them to the observations. Our simulated spectra manage to 
fit both the X-ray continuum and lines at 4 epochs satisfactorily, while our computed 
light curves are in good agreement with additional flux-monitoring data sets. These 
calculations allow us to infer the nature and structure of the circumstellar medium, 
the evolution of the SN shock wave, and the abundances of the ejecta and surrounding 
medium. The data imply that SN 1996cr exploded in a low-density medium before 
interacting with a dense shell of material about 0.03pc away from the progenitor star. 
We speculate that the shell could be due to the interaction of a blue supergiant or 
Wolf-Rayet wind with a previously existing red supergiant (RSG) wind. The shock 
wave has now exited the shell and is expanding in the medium exterior to it, possibly 
the undisturbed continuation of the dense RSG wind. The narrow lines that earned 
SN 1996cr its Iln designation possibly arise from dense, shocked clumps in the GSM. 
Although the possibility for an LBV progenitor for this Type Iln SN cannot be com- 
pletely excluded, it is inconsistent with much of the data. These calculations allow us 
to probe the stellar mass loss in the very last phases (< 10* years) of a massive star's 
life (> 10^ years) , and provide another means to deducing the progenitor of the SN. 

Key words: circumstellar matter; methods: numerical; techniques: spectroscopic; su- 
pernovae: individual: SN 1996cr; stars: winds, outflows; X-rays: individual: SN 1996cr 



1 INTRODUCTION 

Core-collapse supernovae (SNe) arise from stars with zero- 
age main-sequence masses > SM©. However, decades of 
research has failed to determine a direct relationship be- 
tween the various SNe types (as classified based on their 
optical spectra and light curves) and the progenitor stars 
that gave rise to the SNe. It is not clear how the type 
IIP, IIL, lib, Iln, lb and Ic types relate to the properties 
of their respective progenitor stars. The progenitors them- 
selves, as it turns out, are not very well known, with the 
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handful of identified ones appearing to be associ ated al- 
most exclusively with Type IIP SNe (|Smarttll2009l ). It had 
been surmised that the main progenitors of core-collapse 
SNe were red supergiants (RSGs) and Wolf-R ayet (W-R) 
stars (|Falk fc Arnetti[l977l : |Podsiadlowskill 19921 ). The explo- 
sion of SN 1987A revealed t hat blue supergiants (B SGs) 
could also be SN progenitors (^Sonneborn ot al. 1987l) , per- 
haps in a binary system ([Morris & Podsiadlowski 120071 : 
IPodsiadlowski et al.l 120071 ). In the last decade there has 
been disc ussion of LBV stars being the progenitors of Type 
Iln SNe (IChuet al.|[l999l:ISalamancall200alKotakfc VinkI 



20091 : 



2006"; 'Vink' '2008': 'Smith| j2008l : iTrundle et 80112008 
Gal- Yam & Leonard 2009|), which is problematic because 
stellar theorists have mainly placed the LBV stage as an in- 
termediate post-main sequence stage, not as a pre-explosion 
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phase (ISchaller et all Il992l: iLangej Il993l: iLanger et all 2 SN 1996CR 



1994 IStothers fc ChinI Il996l: iGarcia-Segura et all Il996l: 
Maeder et al.ll2005l : iMaeder fc Mevnetll2008l ). 



The problem is clear - we rarely know the progenitor 
star that led to a SN explosion, because it has to be typed 
from pre-explosion images, leading to significant ambiguity 
and potential bias. The expansion of the SN shock wave 
and the resulting emission due to circumstellar interaction 
l|Chevalier fc Frans son 1994) opens up another window into 
the exploration of the pre-SN star. The thermal emission 
from this interaction, including the X-ray and optical emis- 
sion, and to som e extent the non-thermal radio emission 
l|Chevalierlll982bl ) . depends directly on the external density. 
Thus an accurate analysis and interpretation of this emis- 
sion acts as a probe of the density profile. In the case of 
core-collapse SNe, which lose a significant amount of mass 
prior to collapse, the surrounding medium is formed by ma- 
terial from the pre-explosion star. Decoding the structure of 
this circumstellar medium therefore will allow us to probe 
the mass- loss parameters of the pre-SN star, which can then 
be linked to the stellar parameters. Thus this provides a way 
to explore the stellar parameters even after the star ceases to 
exist, and allows us to probe the pre-SN properties of classes 
of SN that have heretofore lacked progenitor counterparts. 

In this paper we pursue this method for the unusual 
Type Iln SN 1996cr. This SN, which exploded around 1996 
but was only discovered about 11 years later, is only the 
second SN after SN 1987A which shows increasing radio 
and X-ray emissio n over a susta i ned p eriod of a few years 
(JBauer et al.ll2008h . [Bauer et all (|2008l ') suggested that the 
increasing X-ray and radio emission may be associated with 
a dense shell of material which the expanding SN shock 
wave interacts with a couple of years after explosion. Herein 
we carry out hydrodynamical simulations to evaluate this 
hypothesis, and to decipher the detailed structure of the 
circumstellar medium into which the shock wave from SN 
1996cr is expanding. We compute the X-ray lightcurves and 
X-ray spectra using non-equilibrium ionization conditions, 
which we compare with high-resolution Chandra observa- 
tions. We achieve a detailed agreement which, given that we 
are using a single model across multiple epochs, affirms the 
validity of our hydrodynamical model, and allows us to con- 
strain a range of abundances for both the material ejected 
in the explosion and the surrounding circumstellar medium. 

The plan of this paper is as follows: In ^ we review 
the observational details of SN 1996cr, including our recent 
High Energy Transmission Grating (HETG) spectra which 
motivated this analysis. Sj3] describes the reasoning behind 
our model of the circumstellar medium (CSM), and the re- 
sults of the hydrodynamical modelling of the SN ejecta in- 
teracting with the CSM. Our techniques for computation of 
the X-ray light curves and spectra from the hydrodynamical 
models, and the resultant X-ray emission, are outlined in 21 
§5 puts the narrow lines that earned SN 1996cr its Type Iln 
designation into the context of the overall model. In fjB] the 
implications for the progenitor, and the abundance determi- 
nations, are discussed in depth. Finally, ^summarizes our 
results and outlines future work in this area. 



This luminous type Iln SN was only identified in the n earby 
Circinus Galaxy (3.7±0.3 Mpc; iKoribalski et al.ll2004r) ~11 
y ears after it was believed to have exploded IJBauen 120071 : 
Bauer et al.ll2008l ). but its remarkable evolution was fortu- 
itously captured in archival data from HST, Chandra, XMM- 
Newton, Spitzer, and several ground-based optical and radio 
observatories taken for other purposes. Initially undetected 
at radio and X-ray wavelengths, the SN emission jumped by 
a factor of at least '^ 300 at radio wavelengths about 2 yrs 
after explosion. It continued to increase modestly for 8 years 
before finally exhibiting signs of a potential X-ray and radio 
turnover in late 2008. 

Fig.[T]shows the radio and X-ray evolution of SN 1996cr. 
Unfortunately, only upper limits exist for the first few years 
after the explosion. But even these are sufficient to illustrate 
the low level of emission in the initial years, followed by a 
steep rise in the radio, and a less steep but no less signifi- 
cant rise in the X-ray emission, which lasted for roughly a 
decade. This behavior is reminiscent of the more fa mous SN 
1987A, wh ose emission follows a similar pattern IjMcCravl 
|2003", 2007), albeit with a factor of ~ 1000 lower luminosity 
(JBauer et al. 200^). Emission from SN 1987A was explained 
bv lChevalier fc Dwarkadas! (|l995l ) as arising from the inter- 
action of the SN shock wave with a dense HII region inte- 
rior to the circumstellar ring seen in the equatorial region. 
Further interaction with the denser ring should result in a 
continuously increa sing light curve in the case of SN 1987A 
(|Dwarkadaj|2007bD . 

Early (2004) HETG X-ray data of SN 1996cr showed 
hints of the complex velocity structure, similar to that seen 
in the optical. This, coupled with the similarity with SN 
1987A and the need to investigate such an unusual SN fur- 
ther, motivated a deep 485ks Chandra HETG observation of 
SN 1996cr in early 2009, with regular monitoring thereafter. 
A companion paper (Bauer et al. 2010, in preparation) de- 
scribes the analysis of the HETG spectrum. In this paper we 
concentrate on the interpretation of the X-ray emission, and 
the computation of the structure of the medium into which 
the SN shock wave is expanding. The steep rise in the X- 
ray and radio, and the more recent turnover in the emission 
(see Fig. [T} point to an interaction of the SN shock wave 
with high density material, from which the shock wave has 
now emerged. This possibility was alre ady indicated by a n 
analysis of the radio spectrum earlier in I Bauer et al.l l|2008r ). 
Herein we provide a more quantitative assessment of the 
various properties of the circumstellar medium based on the 
X-ray data. 

The 2009 HETG spectrum (shown later in Fig. |6} is 
probably one of the most detailed X-ray spectra of a young 
supernova, next only to SN 1987A. The well-resolved emis- 
sion lines allow us to estimate the relative abundances of the 
ejecta to reasonable accuracy. Combined with lower signal- 
to-noise HETG spectra from Chandra in 2000 and 2004, and 
a high signal-to-noise XMM-Newton p-n spectrum in 2001, 
as well as other CCD-quality X-ray data, it provides a strong 
basis for accurate comparison of our simulated spectra with 
the observations, allowing us to fine tune our models. 

Any model of the medium into which SN 1996cr is ex- 
panding must be able to explain the increase in the radio and 
X-ray lightcurves, and the spectral details. The model that 
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Figure 1. (Top) The radio liglit curve of SN 1996cr. All available data from 843 MHz to 89 GHz are shown (some of which is previously 
unpublished and will be presented in Bauer et al., in preparation). Only upper limits are available during the first 3 years. The radio 
flux clearly shows a sharp increase starting from around 700 days, a flattening of the light curve around 2000 days, and a decrease and 
turnover around 5000 days. On top of this luminosity evolution in the 8.5GHz band, we see a gradual decreasing absorption at lower 
frequencies (1.4GHz) that leads to a smoother transition between 1000-3000 days. (Bottom) The X-ray light curve of SN 1996cr over 
time. Unlike the radio, once detected the X-rays show a more-or-less constant increase over time, and no flattening. However the late-time 
turnover is present in the X-ray emission also. The upper limits in the first 3 years are quite high, with the points around 700 and 900 
days being the most constraining. Sample power-law fits to the light curves are also shown. The simulations in this paper are geared 
towards explaining the X-ray emission, while being consistent with the radio behavior. 



we put forth in this paper is constrained by the following 
observational parameters: (1) The range of explosion dates, 
from March 1995 to March 1996 (sec Bauer et al. 200i). (2) 
Observed X-ray upper limits prior to year 2000 and detected 
fluxes thereafter (Figure [1]). (3) A radio VLBI measurement 
that provides a size measurement, which we take to be the 
radius of the outer shock. This gives the radius of the outer 
shock as about 2.8 x 10^^ cm at 12 years (with a statisti- 
cal error of « 20%). (4) Emission- line velocities seen in the 
X-ray and optical spectra, ranging from 2000-5000 km s~^. 
(5) Fits to the temperature (kT) and the absorption column 



density at various times, obtained from the X-ray spectra at 
several epochs. (6) Fits to prominent lines seen in the emis- 
sion spectra, which constrain CSM and ejecta abundances. 
(7) The radio hght curve (Figure [l|. 



3 EJECTA AND CSM MODELS AND 

HYDRODYNAMICAL CALCULATIONS 

In order to study the interaction of the SN ejecta with 
the surrounding circumstellar medium, we need to delineate 
the density profile of the ejecta as well as the surrounding 
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medium. In this first paper we make the assumption that the 
density structure, of both the ejected material and the sur- 
rounding medium, is spherically symmetric, and therefore 
that a one-dimensional (ID) model will suffice. Although 
the complexity of line emission seen, especially in the optical 
spectrum, may indicate otherwise, a spherically symmetric 
model is a reasonable starting point, given the lack of any 
imaging constraints. The validity of such an assumption can 
only be determined by computation of the light curves and 
spectra resulting from these calculations. As we shall show 
later, the spherically symmetric assumption appears to be 
adequate in this regard. The companion paper by Bauer et 
al (2010) will explore the line shapes and emission in detail 
for deviations from spherical symmetry. This work will then 
serve as a launching pad for more complicated explorations 
of the density structure. 

Following the work of IChevalien (|l982al ) and 
IChevalier fc FranssonI ([199J), we denote the ejecta density 
dropping off as a powe r-law with radius (or velocity) . 
iMatzner fc McKed (|l999l ) narrows this power-law to a value 
of around 9 or 10; in the current work we adopt a value 
of 9. In order to conserve mass and energy, the power-law 
density profile cannot extend all the way to the origin, 
therefore below a certain velocity the density is assumed 
to remain constant. This ejecta profile model therefore has 
two free parameters - the normalization of the outer density 
profile and the transition-to-plateau velocity; alternatively, 
these two parameters can be expressed in terms of the 
ejecta energy and mass. 

Note that for early times, until the reverse shock reaches 
the plateau region, only the single parameter of the density 
normalization is involved in the hydrodynamics and there is 
a degeneracy between ejecta mass and energy. In this case 
one can set the energy to a canonical value and use the 
ejecta mass as the free parameter. C hanging the explosion 
energ y scales this mass accordingly fCheval ier fc FranssonI 
Il994h . Our late time modeling of SN 1996cr shows that 
in order to match the light curve, the reverse shock must 
have passed into the plateau region. Fitting the X-ray spec- 
tra constrains the ejecta energy and mass to be 10^^ ergs 
and 4.45 Mq. The latter mass value is very reasonable 
for the total ejected material in a core-collapse explosion 
l|Wooslev et al.ll2002l ). We note that this gives a characteris- 
tic velocity Va = 2 x Ekin/Mejecta = 4740 km s~^ which i s 
suggestive of a stripped envelope SN (IMaurer et al.ll2009h . 
In ^6.1l we will discuss the connection with stripped envelope 
SNe in more detail. 



3.1 CSM Density Profile -Power-Law Model 

Elucidation of the density structure of the CSM required 
considerable trial and error. The simplest assumption of a 
CSM is one with a power- law density profile r"" , where s = 
for a constant density medium and s = 2 for a wind. A 
wind whose mass-loss parameters are changing with time 
can result in other values of s. Therefore we first explored 
the validity of such monotonic-with-radius models. 

The light curve indicates that the X-ray emission be- 
gan increasing linearly about two years after explosion. In 
theory it is possible to have the X-ray emission increasing 
with time even if the ambient density is decreasing with ra- 
dius. This can be shown in the context of the self-similar 



model l|ChevalieiJll982al ) for SN evolution. Detailed results 
for the X-ray emission using the self-similar model, for var- 
ious values of the parameters mentioned below, and inves- 
tigating several different ass umptions for the forw ard and 
reverse shocks, are given in iFransson et al.l l| 19961 ). Herein 
we present a simplified analysis that encapsulates the im- 
portant ideas relevant to this work. 

The X-ray luminosity Lx from a source depends on the 
electron density Tie, the emitting volume V and the cooling 
function A 



n/AV 



(1) 



We assume that the density of the medium into which 
the SN shock wave is expanding goes as r"". The emission 
arises from a thin shell of radius Ar at a mean radius r , 
whose volume V can be expressed as 47rr^ Ar. Note that 
for self-similar evolution, Ar « r, and therefore 1/ oc r^. 
For a SN in the early stages, the temperature is going to be 
much larger than lO^K, and t he cooling function is assumed 
to vary as T^'^ (|Chevalier fc Fransson 1994[ ). For a strong 
shock, T (X Vs^, and therefore the cooling function A (x Vs (x 
r/t in the self-similar case. 

Therefore we get 
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For s=2, a SN expanding in a wind with constant mass- 
loss parameters, this gives the well-known result that the 
emission decreases inversely with time, L^ oc t~^ . 

For s=l, we get L^ oc r^ /t. In the self-similar case, 
r oc t", and therefore we get that L^ oc t^°'~^ . The parameter 



Q = (p -3)/(n-l) for s=l with psiv = At~^v~^ (|Chevaliej 
Il982al ). which is greater than 0.5 for all n > 5. Therefore 
2a — 1 > 0, and the power-law exponent is always positive, 
leading to an X-ray evolution that increases with time, even 
though the density is decreasing as r~^. 

For s=0, Lx oc i''"-\ With a > 0.4, (0.4 being the 
value it would have in the Sedov- Taylor stage), the X-ray 
luminosity increases with time for all values of a. 

Theoretically, it is therefore possible to envision an in- 
creasing X-ray luminosity with a CSM whose density is con- 
stant or even decreasing with radius. The radius can then 
be expressed in terms of the self-similar solution. We use the 
self-similar solution to show that this approximation is not a 
good one in the current situation, for the following reasons: 

(i) The radio spectra indicate t hat the absorption was 
high even two years after explosion (JBauer et al.ll2008r ). and 
that as the SN shock expanded within this region the free- 
free absorption continually decreased. This suggests a region 
of high absorption into which the shock was expanding two 
years after the explosion. The absorption decreased as the 
shock made its way through this region. It would be difficult 
to explain this with a constant or continually decreasing 
density profile. 

(ii) In the discussion above the post-shock temperature 
is related to the shock velocity, which varies as Vs oc t°'~^ . 
Thus the temperature will vary as T oc i^"~^. For example. 
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in order to get an X-ray luminosity L^ oc f^", roughly inter- 
mediate between the hard and soft luminosity slopes seen in 
Figure [l], we can take s — 0.6 and n = 9. Then a — 0.7142, 
and Lx oc f^" from equation [S] This would mean that the 
temperature would vary as t""'^^, which equates to an ex- 
pected temperature drop of 1.7 between 2000 and 2009. 
Such a temperature variation is not seen. Instead the X-ray 
data (Bauer et al. 2010) above 2 keV show little spectral 
variation, and hence little apparent temperature evolution, 
over 9 years. The spectrum below 2 keV does vary somewhat, 
although it is not clear whether this is because of decreasing 
absorption in the soft X-ray regime or whether it is because 
the spectrum is actually becoming softer. 

(iii) A radius of 2.8el7 at a time of 12 years implies an 
average velocity of 7500 km s~^. The optical and X-ray spec- 
tra do not show velocities exceeding 5000 km s~^, and in 
most cases somewhat lower. If we assume spherical symme- 
try and an average velocity of 4000 km s~^ from 2000 to 
2008 (when spectra are available) then it suggests that the 
average velocity in the first 4-5 years (depending on explo- 
sion date) must have exceeded an average velocity of 12,000 
km s~^. A power law model has shock velocity varying as 
Vs OC t"~^ , which would make it extremely difficult, if not 
impossible, to fit such a velocity profile for reasonable val- 
ues of a. Whereas a model with a lower CSM density early 
on and a higher density later could fit it easily. Using our 
earlier example, we can construct a model with n = 9 and 
s = 0.6 using the self-similar solution, which satisfies the 
constraint that the shock radius be 2.8el7 at 12 years. This 
model gives a velocity of 5290 km s~^ at 12 years. While this 
is somewhat high, it could be considered acceptable within 
the error bounds. However such a model would also predict 
a velocity of almost 6000 km s~^ at 8 years and 6450 km 
s"'^ at 6 years, which is not supported by the observations. 
Note that a constant density model with n=9 would have 
Vs OC t'"'^"^, which would make this discrepancy worse. 

(iv) In a model with monotonically decreasing density, 
the density of the CSM would need to be significantly high 
throughout the evolution for the shock velocity to decelerate 
to a few thousand km s~^ in a few years, as suggested by 
the X-ray and optical spectra. For instance, in the above 
model with n = 9 and s — 0.6, the density at 5 x 10^® 
cm would be 7.71e-21 g cm~"^. This is about 3800 times 
greater than what is in our simulations presented in the next 
section, which would result in the X-ray emission, especially 
the hard X-rays, being about 10^ times greater on average, 
significantly exceeding the upper limits derived from archival 
observations. 

(v) The radio emission is increasing significantly with 
time after about 2 years (Figure [l]), which is hard to ac- 
complish in a model with decreasing density. In particular 
free-free absorption by itself cannot account for this increase 
l|Bauer et al. 2008), there still needs to be a flux increase 
> 100. This can best be implemented with a density in- 
crease, as the radio emission, althoug h non-thermal, indi- 
rectly depends on the external density (|Chevalieilll996l ). 

(vi) The X-ray emission increases significantly after two 
years. A simple power-law increase in the X-ray flux, as pre- 
dicted by the self-similar model, cannot match the data, 
which requires at least two distinct power-laws, with the up- 
ward slope becoming much steeper after 2 years. This would 
then require at least two different density slopes, with the 



density decreasing slower in the second case. When com- 
bined with the VLBI radius constraint, and the velocity and 
temperature constraints outlined above, it would be almost 
impossible to make the power-law model consistent with all 
the available data. 

The above arguments suggest that a CSM density pro- 
file that shows a simple power-law decrease would not work. 
Although a self-similar solution is not applicable, many of 
the same arguments can be used to show that an increasing 
power-law density profile from the origin (the stellar sur- 
face) also would not work. Instead, a density profile where 
the density changes after 2 years is suggested. 



3.2 CSM Density Profile - High Density Shell 

The above discussion clearly argues for these basic require- 
ments of any CSM model for SN 1996cr: a CSM whose den- 
sity is low close to the star, then increases sharply with ra- 
dius. This would fulfill both the velocity and density cri- 
teria outlined above, as well as lead naturally to the ob- 
served radio spectra and light curves. One way to arrive 
at such a distribution is with a low density wind medium 
followed by a thin, dense shell, as is often seen around 
mass i ve stars (IWeaver et al.l 1 19771 : IChul l2003l : ICappa et al] 



I2OO3I: IChul bOOSl). and in the nebula around SN 1987A 



( Blondin fc Lundgvisjl 19931 ). In this paper we consider such 
a distribution. The inner and outer radius of this shell, the 
shell density, and refinement of the CSM density structure, 
are obtained via an iterative procedure that consisted of 
adopting a value for the shell parameters, computing the 
hydrodynamic interaction, calculating the X-ray light curves 
and spectra, and re-iterating until a reasonable fit was ob- 
tained between the simulated and observed spectra and 
light-curves. This took on the order of a dozen iterations 
before convergence was obtained. 

The model that we finally converged on has a shell 
with density around 1.28 xlO~^^ g cm~"^ extending from 1 
xlO^^ to 1.5 xlO^^ cm. The X-ray light curves and spectra 
are highly sensitive to the shell parameters, i.e. the inner 
and outer shell radii and shell density, and so manage to 
constrain them quite tightly. Even a 20% deviation of the 
shell radii from these values destroys the agreement with 
the lightcurve, while a factor of 1.3 in the density has about 
the same effect. Therefore the shell itself is strongly con- 
strained by the X-ray data. A rearrangement of the mass 
in the shell, i.e. a modest change in the density profile that 
does not alter the total mass (Figurc[2)), can result in modest 
changes to the early-time X-ray flux, as shown in Figure [5] 
For the CSM inside of the shell, we have assumed a structure 
resemblant of a wind-blown bubble, with an initial freely ex- 
panding wind whose density drops as r~^ ending in a wind 
termination shock, which is assumed to be a strong shock. 
Beyond this is a constant density shocked wind region. The 
wind density, and location of the wind termination shock, 
are constrained only by the density of the shocked wind re- 
gion, and are obtained as described in the Appendix. For 
the standard case described here and shown in Figure (2] the 
value of the wind parameter M/v^ is 6.76 xlO^^g cm~^, 
the wind termination shock is located at 1.63 xlO^® cm, 
and the shocked wind region has a constant density of 8.15 
xlO"^^gcm"^. 
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The structure of the medium beyond the dense sheU 
is difficult to estimate properly, since the shock front has 
not had much time to expand within it. One assumption 
is that the dense shell consists of fully swept-up wind ma- 
terial. In this case the CSM density structure resembles a 
wind-blown bubble formed by the interaction of a super- 
sonic wind with the slower wind from a previous epoch; the 
dense shell is composed entirely of the slower wind mate- 
rial that has been swept-up. The mass of the shell must 
equal the total wind mass up to the outer radius of the shell, 
or Msheii = (M^/Vn,) Rsheii, wherc the subscript w refers 
to the outer wind. Since the shell outer radius is known, 
and the mass is defined once we have the shell boundaries 
and shell density, this provides us with the wind parame- 
ters Mw/Vw Assuming constant wind parameters, the den- 
sity of the material outside the shell at any given radius 
r can then be expressed as that of a freely expanding wind 
p = Mw/Vw X l/(47rr^). This was one of the two assumptions 
that we have tried, and since it gives satisfactory results over 
the short period of time that the shock is interacting with 
this region, we have continued to use it. (The other assump- 
tion involved a constant density profile exterior to the shell. 
However in that case the shell mass was always found to 
be larger than the mass of the swept-up material given that 
density, suggesting a complicated origin for the shell. A con- 
stant density profile either gives too small a density close 
in or too large a density further out, depending on its mag- 
nitude. Furthermore, the total absorption column from any 
constant density profile that provides adequate X-ray emis- 
sion beyond the dense shell increases too fast, and thus such 
a profile cannot be sustained beyond 0.4pc in any reasonable 
model. Since the wind profile anyway fits better, we did not 
explore the constant density assumption further). 

The final density profile for the SN ejecta and CSM 
that we have used, which forms the initial conditions for the 
hydrodynamical simulations, is shown in Fig. [5] The thin 
solid lines superimposed on the constant density shell are 
intended to show how rearranging the mass to perturb the 
density in this fashion can lead to a modest change in the 
X-ray lightcurve (the dashed gray line in Figure [5]). 



3.3 Hydrodynamic Simulations 

We carried out hydrodynamical simulations to compute 
the interaction of the SN ejecta with the CSM. The sim- 
ulations were carried out using VH-1, a three-dimensional 
finite-difference hyd rodynamic code based on th e Piecewise 
Parabolic Method of lColeUa fc Woodward! |l984 ). Ahhough 
cooling is included in these simulations via the inclusion of a 
cooling function, it was not found to play an important role 
in the ID calculations. However the cooling time of the shock 
within the shell was of the same order of magnitude as the 
evolution time, and it is possible that in multi-dimensional 
calculations, where the growth of instabilities can give rise 
to dense clumps, cooling could prove to be more effectivqj. 



The simulation was started with the initial conditions de- 
scribed above, and the SN ejecta allowed to expand into the 
surrounding medium. 

The evolution of the various dynamical and kinematic 
quantities over time are shown in Fig. [S] The SN ejecta ex- 
pand quickly in the low density medium (see Fig. 2} until 
they reach the higher density shell in about 20 months. The 
shock, which was expanding very rapidly up to this time, 
is instantaneously decelerated by the high density, and its 
velocity drops to about 2000 km s~^ due to the high pres- 
sure. As seen in Fig. [S] the high pressure region between the 
forward and reverse shocks gets highly compressed. The col- 
lision of the ejecta with the high density shell sends a lower 
velocity transmitted shock into the dense shell, and a re- 
flected shock is formed that travels back into the ejecta. The 
reflected shock quickly overtakes the original reverse shock, 
thus thermalizing the ejecta faster than the SN shock itself 
would have done. The transmitted shock gradually increases 
in velocity as the highly over-pressured region slowly depres- 
surizes, and then attains a more or less uniform velocity of 
around 5500 km s"'^. 

The interaction of the shock wave with the high density 
shell lasts for about 4.6 years (from t ~ 1.5 to t ~ 6.1 years.). 
The structure of the interaction region resembles the interac- 
tion of a power-law profile with a constant density medium, 
as expected, and the density decreases from the reflected 
shock to the contact discontinuity. Since the pressure does 
not vary much, the temperature is almost inversely propor- 
tional to the density, and increases outwards from the re- 
flected shock, until it reaches a large value at the contact dis- 
continuity. We note that the structure is comparable to that 
described by the self-similar solut ion for a steep ejecta den - 
sity profile colliding with a wall IjChevalier fc Liang 11983 ). 
The comparison is close but not exact, since the shock in 
this case imparts some motion to t he shell, which is not in- 
cluded in lChevalier fc Liangj (|l989l ). It is also interesting to 
note that, given the presence of the dense shell so close to 
the star, the density into which the forward shock is prop- 
agating, while in the shell, is much higher than the density 
into which the reverse shock is propagating. This is contrary 
to a general expansion in a power-law or constant density 
medium, and results in the temperature behind the forward 
shock (~ 10^) being initially lower than the temperature be- 
hind the reverse shock {^ 10^). Therefore we expect that the 
major contribution to the X-ray emission in the 0.5-10 KeV 
wavebands at this early time arises from the shocked circum- 
stellar material. As the shock continues to expand, the tem- 
perature behind the forward shock increases, whereas that 
behind the reverse shock decreases as the shock expands into 
the increasing steep power-law ejecta density profile. When 
the transmitted shock exits the shell and begins to reform 
in the freely expanding wind around it, the reflected shock 
is then interacting with the constant density ejecta, whose 
density is higher than that of the wind medium. The den- 
sity behind the refiected shock is higher than that behind 
the forward shock, its temperature lower, and the emission 



^ Our X-ray spectral calculation, as well as the formation of nar- 
row optical lines, does in fact require the presence of dense clumps, 
but they comprise at most 1% of the mass. In these scenarios 
cooling is important, and in fact the shock within the clumps is 
assumed to be radiative in order to account for the narrow-line 



emission. However this would affect only a small portion of the 
material. In any case, these clumps are not included in the spher- 
ically symmetric numerical hydrodynamics simulations, and they 
would not affect the general hydrodynamical evolution described 
herein. 
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Figure 2. The initial density profile of the SN ejecta and CSM, used as initial conditions in the hydrodynamic simulations. Starting 
from left and moving outwards in radius, the density profile shows the constant density plateau region of the SN ejecta, followed by the 
steeply decreasing ejecta profile which goes as r~^. Beyond this is the circumstellar medium. The ejecta have already crossed the freely 
expanding wind region (not shown). Going outwards in radius, we encounter a uniform density shocked-wind region, the high density 
shell, and the freely expanding outer wind. A suitable temperature was assigned to each component. The pressure (and temperature) 
that we use for the initial conditions is not important, since it is significantly smaller than the pressure of the gas shocked by the SN 
blast wave. The labels attributed to the various sections of the CSM reflect the discussion in ij6] The thin solid is intended to show how 
a rearrangement of the mass in the shell in this fashion, without changing any other parameters, can lead to a modest increase in the 
light curve at early times (< 3 years), shown as the dotted line in Figure [5] 



becomes dominated by the reverse shocked ejecta. It was 
found that a good fit to the late-time emission required that 
the reverse shock enter the plateau ejecta region. The lo- 
cation of the plateau, the velocity below which the ejecta 
density remains constant, was adjusted to match the obser- 
vations, fixing the total ejected mass to the value given in 
the previous section. 

Once the transmitted shock exits the shell, it is expand- 
ing in the lower density wind, whose density decreases as 
r~^. The radius of the forward shock at ~ 12 years is consis- 
tent with that obtained from VLBI imaging to within 5%. 
Thus this approximation has proved suitable so far, but only 
time and future data can confirm whether this approxima- 
tion is valid for the further expansion of the shock wave. 
While most of the X-ray emission at this point is expected 
to arrive from the reverse shocked ejecta, the radius and 
velocity of the forward shock, and therefore the size of the 
remnant, are determined by the interaction with this wind. 

It is interesting to note that although the shock is inter- 
acting with an r~^ wind, the shock structure is resemblant 
of the shock interacting with a constant density medium. 
This is to be expected as it will take a few dynamical times 



for the shock to r each the self-similar solution in the wind 
(jDwarkadasI l2005l l . Over time therefore the entire density 
structure of the shocked interaction region is predicted to 
change, with the density at the contact discontinuity going 
from relatively small to relatively large values. Such a change 
should be reflected in the X-ray emission over the next few 
decades. 



4 COMPUTATION OF THE X-RAY EMISSION 
AND COMPARISON WITH OBSERVATIONS 

4.1 Non-equilibrium Ionization and Lagrangian 

Shells 

Our goal is to compute the X-ray emission from the output 
of the hydrodynamic simulation efficiently and with rea- 
sonable accuracy. Plasma that has been recently shocked 
and whose densit y is low is presumab l y not in i oniza- 
tion equilibrium (IHamilton et al.l 119831 : iBorkowskil |2000| : 
iBorkowski et al.ll2001^ : as a rule of thumb this is the case 
if r = Uet < 10^^ s cm~^ where n^ is the electron density 
and t is the time since the plasma was shocked. Since the 
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Figure 3. These panels depict the hydrodynamic evolution of various quantities: ejecta density (black); CSM density (blue); temperature 
(red); and fluid velocity (green). The grey vertical lines denote the borders of the mass-shells into which the domain is divided up, for 
the purpose of computing and tracking the fluid history. The thick grey line shows the initial density profile to enable the reader to see 
how it changes with time. The last four panels approximate the times when the spectra outlined in §4 were taken. The various panels 
show, from left to right and top to bottom (1) the initial interaction of the SN ejecta with the shocked wind, and the formation of a 
forward and reverse shock (2) The interaction with the dense shell, the compression of the interaction region between the shocks, and the 
formation of a resultant transmitted and reflected shock (3) (HETG-00) The transmitted shock is advancing slowly through the shell, 
while the reflected shock goes up the steep ejecta incline. Note that the interaction region is slowly growing bigger. (4) (XMM-01) The 
forward shock has reached the edge of the shell, and will soon start to expand out into the external wind. (5) (HETG-04) The reflected 
shock reaches the constant density plateau region of the ejecta. Beyond this point, the solution is no longer self-similar, if it indeed ever 
was given the various transitions. (6) (HETG-09) The evolution of the forward shock in the external wind, while the reverse shock is 
expanding back into the ejecta plateau region. 



plasma impacted by the SN shock wave has a density typi- 
cally less than 10^ cm~'^, the time to reach ionization equi- 
librium is measured in years and the X-ray emission must 
be computed using non-equilibrium ionization (NEI) condi- 
tions. In general this is done by calculating the evolution of 
ion fractions using ionization and recombinatio n rates under 
appropriate simplifyin g assumptions (Hamilto n et al.lll983l : 
iBorkovyski et al.lll994 ). Here, however, our approach is to 
make use of existing NEI codes and appropriately interface 
them to the hydrodynamics output. 

We use NEI emission models available in the XSPEC 
package, vgnei and vpsh ock, along with "N EI version 2.0" 
atomic data from ATOMDB l|Smith et al.ll200ll) which has been 
updated to include relevant inner-shell processe^j. These 

^ Files provided by K. Borkowski and online at 
|http : //space . mlt . edu/home/dd/Borko«ski/ | 



codes encode the NEI state with only a few parameters; 
we have chosen to sacrifice some accuracy for speed and 
simplicity. To calculate the X-ray emission from an indi- 
vidual fluid element these models need information that is 
based on the element's history, in particular the ionization 
age, r = Jne{t)dt, the current electron temperature, Te, 
and the ionizatio n-age-averaged temper ature, < kTe > — 
jTe{t)ne{t)dt/T JBorkowski et al.ll200ll ). These quantities 
can be calculated from the set of simulation time steps pro- 
vided the location of a given fluid element can be determined 
and followed in time. To this end we regroup the Eulerian 
output at each time step in a Lagrangian manner, using cu- 
mulative mass as the Lagrangian parameter to deflne and 
track a set of mass intervals, representing "shells" in 3-D. 
These mass shells, of order 50 of them, are not as finely 
spaced as the hydro radial bins (thousands), so we calculate 
shell-averaged hydrodynamic fluid parameters: pressure pj , 
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Figure 4. The velocity of the SN forward shock during the evolution. It is extremely high (exceeding 20,0000 km/s) in the first few 
months, but then drops dramatically when the shock encounters the higher density shell, and the interaction region is highly compressed. 
As the compressed region slowly expands, the velocity rises again as the built-up pressure is released (A similar effect was observed in 
SN 1987A). After about 6 years, the shock exits the shell, and the velocity increases as it encounters a lower-density medium. As the 
density profile of the shocked interaction region changes from that in a constant-density medium to that in a wind medium, the velocity 
remains more-or-less constant, at a velocity around 5500 km s~^. 



density pj , and radial velocity Vj , for each shell at the time 
steps, tj. 



4.2 Abundances, fi, and Te 

The hydrodynamic variables are converted to values rele- 
vant for the X-ray emission, rie, Te, through the abundances 
and the related mean mass per particle, jx. We assign an 
abundance to each shell of material: shells interior to the 
contact discontinuity are assigned ejecta abundances and 
those exterior are CSM. It is possible (and as we argue in 
the discussion section, most likely) that the CSM can be di- 
vided into two phases, the portion interior to the shell and 
the portion including the shell and exterior to it. However 
the evolution of the SN ejecta into the region interior to the 
shell is not constrained by the data. Therefore a single set 
of abundances is assumed for the CSM for simplicity. It is 
likely that the ejecta may be layered, as is expected for a 
massive star, and different abundance distributions should 
be expected for different layers. However, we have no a priori 
indication of any of the abundance distribution, and 4 spec- 
tra of different resolution spread over 9 years, only the last 
two of which are actually tracking reverse-shocked emission 
according to our calculations, would be unable to distin- 
guish between abundances of various layers. Therefore as a 
first approximation we have considered a single set of ejecta 
abundances, although our code is easily capable of assigning 
a unique abundance to each shell of material if required. 
Define X{Z) as the relative abundance values with re- 



spect to a refere nce number-ratio abundanc e set, Aags9(Z); 
here we use the lAnders fc Grevessd (|l989l ) abundance set. 
The mean atomic weight can be calculated as: 



MA 



J2zmzX{Z)AAG89iZ) 

Y,^X{Z)Aags'AZ) 



(4) 



where mz is the mass of an ion of element Z. Similarly the 
mean charge state averaged over all ions is given by: 



lA = 



J2zQziTcm)X{Z)AAG8u{Z) 

J2^X{Z)AAGS9iZ) 



(5) 



where Qz{Tcie) gives the average charge state of ion Z, i.e., 
the average number of free electrons per Z-ion, and Tcie 
refers to the temperature in collisional ionization equilib- 
rium. This value changes with progressive ionization, how- 
ever it is dominated by the low-Z elements which are quickly 
ionized, so we approximate its value based on ionization 
equilibrium at a fixed temperature, Tcie = lO'^ K. 
The mean mass per particle is then: 



M = 



^J.A 



where fi here has units of mass.U The densities of electrons 
and ions can now be obtained from the fluid density as: 



•^ Note that it is common to informally express fi in implied 
atomic mass units, /i/rriamu, e.g., "/i = 0.5 for a fully ionized 
H plasma". 
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na{Z) Aags9{Z) I Y.^ X{Z) AAGS9iZ). 

We can now compute tfie average particle temperature 
in an element of plasma from the hydrodynamic pressure 
and density: fcThydro ~ l-i-p/p- Since the shocks are coUi- 
sionless, most of the post-shock energy is transferred to 
ions, and the electron temp erature is general l y a f raction 
13 of the ion temperature (|Ghavamian et alj 120071 ). This 
is accounted for by setting Te = /9ah('fsh ) Thydro with the 
Bsh_vahie set to be « vj^ as given in iGhavamian et alJ 
l|2007l ): the shock velocity is calculated self-consistently from 
fcTiiydro = (3/16)/x Vsh^. Other choices for /3sh(^sh) can be 
easily implemented if desired. For example, we clipped /3 to 
have a minimum value of 0.05 as an approximation to the re- 
sults of Ivan Adelsberg et al.l (|2008l '): however, this produced 
no significant changes in our calculated flux or spectra as 
the P values of the mass-cuts are generally well above 0.1. 

Finally, the /3 value is modified to include its evolution 
due to Coulomb interactions: Te/Thydro approaches 1.0 as 
the ionization age approaches requii(i'sii)- The dependence on 
shock velocity is given roughly by: Tequn ~ 10^" (i'ah/i'40o)^, 
with U400 ~ 4x10^ cms~^; this follows from the T'^" varia- 
tion of the Coulomb equilibration time constant. For smaller 
values of r, /3 chan ges such that / 3 ^ (T /Teauji)"'^^, e.g., as 



ifi 



seen in Figure 4 of iMichael et al.l (|2002l '). With these rela- 
tions the value of /3 is increased from its initial post-shock 
value when the r of the mass-cut increases sufficiently. 

4.3 Implementation and Embellishments 



Using the above equations, for each shell at each time 
step we can calculate the input v alues needed for the 
vgnei model (|Borkowski et al.l [2OOII '): kT, Tau, <kT>, and 
the norm= 10"" AV UeUH /(47rd^), with AV the shell vol- 
ume and d the source distance. The vgnei abundance val- 
ues are the appropriate set (ejecta, CSM) scaled so that the 
H abundance is 1. The calculations, from reading the hy- 
dro output files to evaluating the XSPEC spectral model, 
are all done with custom routines written in S-Lang, the 
interactive scripting lan guage used in the ISIS package 
JHouck fc Denicolall2000l ). 

The flexibility of the code also allows us to implement 
additional modifications to improve the fidelity of the emis- 
sion calculation as needed. For most mass shells we use the 
vgnei model with an appropriate r value based on the shell 
history as described above. For the case where a mass shell 
is partially shocked, the shock front is within the shell and a 
constant r value is not appropriate since there is a range of 
r starting at r = up to a maximum value. In this case we 
can calculate the shell's emission using the vpshock model, 
which allows us to set a range of r from a low (zero) to a 
high value. Inclusion of these partially-shocked shells adds 
some line emission of lower ion states and smoothes the light 
curves. 

There are indications that the X-ray emission from SN 
1996cr is not simply the sum of optically thin components: 
the as-fit Nh decreases with observation epoch and the X- 
ray line shapes at high-resolution are asymmetric. These 
data are examined in more detail in the companion paper 
(Bauer et al. 2010, in preparation) but the methods to in- 
clude them in the emission model are briefly summarized 
here. To self-consistently model the Nh effect, we can add 



"internal absorption" proportional to the amount of un- 
shocked CSM material external to the forward shock. In 
the case of the asymmetric line shapes it seems likely, and 
our hydro model predicts, that there is substantial opac- 
ity through the SN ejecta core. Hence some of the red- 
side emission is obscured producing asymmetric lines and 
reducing the observed flux. Based on simple 3D modeling 
(jDewev fc Nob le 2009), a value of 0.75 has been used to ac- 
count for the amount of flux typically obscured by the core, 
and a blue-shift and Doppler Gaussian broadening, propor- 
tional to the hydro shell velocity, are applied as a first ap- 
proximation to the line shape from the shell. 

4.4 Comparison with the Data 

As the previous sections imply, the X-ray emission from a 
given hydro simulation depends on some additional param- 
eters and assumptions. Perhaps the most influential of these 
are the elemental abundances assigned to the hydro plas- 
mas. These abundances show up in two main ways: i) they 
determine the /i value of the plasma and hence affect the 
electron temperature and the number of electrons per ion 
(section 14. 2|) . and ii) the abundances of the higher Z ele- 
ments are responsible for the X-ray emission lines in the 
spectra. Because we use the XSPEC vgnei (vpshock) rou- 
tine, the abundances of the 13 elements currently treated by 
that model are the ones relevant here: H, He, C, N, O, Ne, 
Mg, Si, S, Ar, Ca, Fe, Ni. These abund ances are specified 
relative to the I Anders fc Grevessd (| 19891 ) values. 

The amount of the lowest-Z elements, especially the ra- 
tio of H to He, is most responsible for determining the mean 
mass per particle p and thus changing the electron temper- 
ature proportionally. Temperature changes lead to changes 
in the continuum shape and in the ratios of He-like to H- 
like line emission from higher-Z elements. Because their line 
emission is not clearly visible here, yet they contribute most 
of the X-ray continuum, we refer to the elements H, He, C, 
N, and O as "continuum elements" . Using only X-ray data it 
is hard to decide among different possibilities for these "con- 
tinuum" abundances: H and He produce no lines in the X-ray 
range, and depending on the temperatures and amount of 
absorption involved, it can also be difficult to see lines to 
constrain C, N and even O abundances with X-ray spectra. 
For this reason we have little handle on the abundances of 
C, N, and O in SN 1996cr. 

With the continuum abundances set, the abundances 
of the higher-Z elements can generally be varied quite a bit 
without effecting the /i value and in this sense the high-Z 
abundances are relatively uncoupled to the hydrodynam- 
ics. However, because line emission from these elements, 
Ne through Fe, can be clearly seen in our X-ray spectra 
these abundances are well constrained by the data. Two 
other "degrees of freedom" in generating the observed X- 
ray emission from the hydrodynamics are absorption, both 
foregr ound and internal, and clu mping in the ejecta and/or 
wind IChugai fc Danzigeilll994l 'l. Both of these have little 
effect on the higher energy flux, 2-8 keV, but can produce 
20% to 50% changes in the low-energy flux. In this low 
energy range there is some degeneracy between the CNO 
abundances, the Nh, the clumping parameters, and the Fe 
abundance through Fe-L lines. For our multi-epoch model 
we have set the Nh at the HETG-00 and HETG-09 epochs 
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to 9.6 xlO^^ cm"^ and 7.8 xlO^^ cm~^, respectively, up 
to 6x10^^ cn\~^ of which is Ukely from absorption in the 
Milky Way and Circinus Galaxy disks. We then interpolate 
to other epochs based on the amount of unshocked GSM. 
For the swept-up shell clumping, which likely extends over 
a range of values, we have chosen a density factor of 7.0 and 
a filling fraction of 0.0020; these are sufHcient to generate 
the Si He-like lines seen at early times. 

The output of our calculations provides the emitted X- 
ray fiux versus time in the 0.5-2 keV band and the 2-8 keV 
band that we have used, as well as spectra and spectral 
model parameter files at each epoch. Furthermore, we have 
performed a detailed comparison of the line-shapes between 
the data and model Si and Fe lines at the HETG-09 epoch, 
which are described in the companion paper by Bauer et 
al. (2010). 

Figure [5] shows the comparison between the fiux cal- 
culated from our hydrodynamic simulations and the obser- 
vations. The abundance list used to generate these plots is 
given in Table 1. The separate contribution of the shocked 
ejecta and the shocked GSM to the flux are also shown. The 
modelled flux is within 20% except for the year 2000 data 
where it is low by 30-50 %. As Figure [3] (panel 3, 5.02 years) 
shows, at this time the flux is determined by the inner-half 
of the dense shell. We have verifled (Figure [S]) that shifting 
10% of the dense shell mass from the outer half to the inner 
half can produce the needed fiux increase with little change 
at later times. While it is tempting to expand the shape of 
the dense shell into a set of basis functions and adjust them 
for the best light curve we feel this would be "over tuning" 
our simple model and limited data. Gertainly in future when 
IXO will obtain a deep, well-sampled dataset of a similar SN, 
such an approach could be very fruitful. 

It is clear that for about the first 7 years the shocked 
CSM material dominates the X-ray lines and the continuum, 
but after that the contribution of the ejecta increases as the 
forward shock exits the shell. By about 8 years the shocked 
ejecta becomes the primary contributor. If our approxima- 
tion of a wind medium exterior to the dense shell is correct, 
then the shocked ejecta will remain the primary contributor 
to the total flux for more than 25 years. It is interesting to 
note that this is exactly the opposite scenario as deduced for 
another SN with a high X- ray fiux, SN 1993J JNvmark et all 
l2009l : lGhandra et al.ll2009l ). where the reverse-shocked ejecta 
was found to dominate the X-ray emission for the first sev- 
eral decades. The dominant early contribution of the CSM 
is clearly due to the presence of the bubble and dense shell 
in this case. 

Our method enables us to compute simulated spec- 
tra using non-equilibrium ionization, which we can com- 
pare to the observations. Observed Chandra HETG spectra 
are available for 2000, 2004, and 2009, and observed XMM 
spectra for 2001. One important point of our method is 
that the same hydrodynamieal simulation can produce spec- 
tra that match the observed ones at 4 different epochs, with- 
out any unnecessary free parameters. Given the transition 
from GSM to ejecta dominance in the emission, the GSM 
higher-Z abundances were obtained by using fits to both 
the HETG 2000 and XMM 2001 spectra, while the ejecta 
abundances were set by fitting to our 2009 spectrum. 

Figure [6] shows the comparison of the simulated and ob- 
served spectra at all the 4 epochs at which observed spectra 



are available. Our simulated spectra are able to fit the con- 
tinuum emission very well, as well as most of the major lines 
(H and He-like transitions of the various lines are shown on 
the plot to aid the reader). We also manage to fit the line 
transitions in the hard X-ray band adequately. Initially, our 
simulated spectra did not adequately reproduce the He-like 
transitions of many elements, especially those that occur in 
the soft X-ray band below 2 keV. In order to fit these bet- 
ter, we assumed that there existed a component, about 1% 
by mass, with a density that was about an order of mag- 
nitude higher than that of the surrounding medium, and 
whose temperature was consequently an order of magnitude 
lower in order to preserve pressure equilibrium. High-density 
clumps may be responsible for this component. The rich and 
complicated line structure, at both the X-ray and optical 
wavelengths, is a good indicator of clumping, and the small 
percentage required to fit the various lines is well within the 
realm of possibility. 



5 SN1996CR AS A TYPE UN - CLUMPS AND 
NARROW-LINE FORMATION 

SN 1996cr was classified as a Type Iln supernova, based on 
the presence of nar row lines in the VLT spectrum from 2006 
(JBauer et al.ll2008l ). In the model presented herein, by 2006 
the shock had exited the shell and was expanding in the 
external wind region. Thus the narrow optical lines could 
not arise due to the shock within the shell, whose maximum 
velocity in any case was much higher than the ~ 670 kms~^ 
FWHM of the narrow component JBauer et al.lbOOSi '). It is 
possible that the actual FWHM could be even smaller and 
the lines even narrower, beyond the limit of resolution of 
that particular observation. 

Gould the narrow lines arise from recombination in the 
unshocked, ionized GSM, which we suggest (below) is a RSG 
wind? The velocity of the putative RSG wind would be low, 
< 50kms~^, and this line would not be resolved. In order 
to get the observed fiux in the narrow lines, the RSG wind 
must have a certain m inimum mass. Using equation 11 in 
ISalamanca et al.l (|l998l ). we find that in order to get the 
required H/3 luminosity, a H number density in excess of 
10* cm~"^ is required at the position of the shock, at the 
time the VLT spectrum was taken (Jan 2006). This is more 
than 3 orders of magnitude higher than the density in our 
simulations, suggesting that in our scenario this line could 
not have come from the unshocked GSM. 

It is possible that the lines arise from shocked dense 
clumps suggested in i]4.4l In order to see if this explana- 
tion is feasible, and to estimate the density, number and 
filling fraction of clumps, we follow the method outlined by 
IChugai fc Danzigeil ([199J). From the simulations, the aver- 
age pressure in the shocked region at 10.5 years is about 3.9 
xlO~^ dynes cm~^. If the shock transmitted into the clump 
has a maximum velocity Vc of 67 km s~^ (the appr oximate 
velocity of the narrow Ha line; iBauer et al.l 120081 ) and is 
in pressure equilibrium, then the clump density pc is 8.69 
xl0~^^ g cm~^. Using the average value of the mean molec- 
ular weight for the GSM (see Table 1), this would mean a 
number density Uc ~ 6.9 x lO^cm"^. This density is consis- 
tent with the "clumps" assumed in the previous section. 

The cooling time of the shock in the clumps must 
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Figure 5. The X-ray flux in tiie soft (0.5-2 keV; blue lines) and hard (2-8 keV; red lines) bands computed from our simulations, and 
compared to the observations (circles - error bars are shown). The data are as observed, not corrected for Njj absorption. The emission 
from the higher density "clumps" is included, although it does not make a significant contribution to the light curve, especially at late 
times. An explosion date of 1995.4 is used. The individual contributions of the shocked ejecta (dotted lines) and shocked CSM (colored 
dashed lines) to the hard and soft fluxes is also shown. The CSM dominates the early evolution up to about 7 years, then the shocked 
ejecta begin to take over and dominate the later evolution. The abundances used for the figure are given in Table 1. The long-dashed 
grey lines indicate the light curve that would be obtained if the density profile depicted by the thin solid line in Figure [2] was used, 
showing how a slight rearrangement of the mass in the shell can affect the early-time light curve. 



be small compared to the destruction time for the clumps 
by Rayleigh- Taylor and Kelvin-Helmholtz instabilities. The 
tem perature of the sho c k driv en into the clumps, follow- 
ing [cMgiL&iDaniigei l|l994 ). is Tc = 1.36 X 10^ vl^ ~ 
6.1 X 10^ K. The cooli n g func tion at that temperature, from 



IChevalier fc FranssonI (| 19941 1. has a value of 5.26 xlO" 
ergs cm'' s~^. Assuming that the cooling time is less than 
the time-scale on which the instabilities develop (~ ac/vc) 
gives a minimum size for the clumps of ac ~ 1.55 x 10^^ cm. 
As suggested by IChugai fc Danzigeij ([199J) , the nar- 
row component arises in the cool material behind the clump 
shock, which radiates due to the reprocessing of X-ray emis- 
sion from the hot shocked gas. As the shock cools, it will 
progressively radiate X-ray, then UV and then optical emis- 
sion from the cooling region behind the cloud shock. An 
upper limit to the luminosity of the ensemble of Nc shocked 
clumps is given by: 



0.25 n al pc v^ N^ 



(7) 



The luminosity of the narrow component of Ha is 2.79 
Using the values of pc and Vc above gives 



X 10^* ergs s 



iVc~ 1359 0^^5 (8) 

where Ucis is the clump size in units of 10^''' cm. Keeping 
in mind the lower limit on ac, we adopt ac ^ 10^" cm, 
which gives Nc ~ 13.6. The total mass in clumps is about 



0.025 A/0. Note that the (3D) filling fraction of the clumps, 
fc ~ 1.2 X 10~^, is quite low, and consistent with that ob- 
tained from the X-rays given the large uncertainties. The 
(2D) covering factor ~ 0.027, thus indicating that they will 
not obscure the SN. However these clumps could add addi- 
tional column density along certain sightlines. 

A consistency check on the above calculations can be 
done by comparing the ratio of the velocities. Since the 
clumps are in pressure equilibrium with the interclump 
medium, the ratio of the interclump to clump shock veloc- 
ity must be proportional to the square-root of the clump 
to interclump density. The density of the unshocked gas at 
this epoch in the simulation is 1.34 xlO"^" g cm~^. Given 
the density of the clumps mentioned above, the jump in 
density between the clumped and interclump medium is a 
factor of 65. If these are in pressure equilibrium, then the 
velocity ratios must be equal to the square-root of this, or 
about 8.06. Therefore, if the velocity of the shock within the 
clumps is 670 km s~^, this implies that the shock velocity 
must be about 5400 km s~^. This is close to the velocity of 
the CSM shock in our simulations (see figure Q, as would 
be expected. However, this is not the vel ocity of the broad 
Ha line mentioned in iBauer et al.l (|2008l ). which would be 
assumed to reflect the shock velocity. The inference is that 
either the broad Ha velocity has been underestimated, or 
that it does not arise from the shocked CSM, but perhaps 
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Figure 6. Comparison of tlie observed (black) and calculated (red) spectra at various epochs. H and He-like transitions of various 
elements are marked on the plot. Our simulated spectra arc able to match the continuum beautifully, and the comparison with various 
lines is adequate. A single normalization factor is adjusted at each epoch - the spectra do not contain any other free parameters. Further 
details of how the spectra were calculated are given in §4. 



from the shocked ejecta. Computing the FWHM of Hq de- 
pends on fitting the continuum precisely, which depends on 
the degree and nature of the reddening towards the SN, 
which is not accurately known. It is possible that the Ha 
has been underestimated, and that a more accurate fit to 
the Ha line will yield a FWHM compatible with this value. 

The computed values for the number, size and 
density of clumps are c omparable to those found by 
IChugai fc Danzigeij (|l994l ). and are consistent with the 
"clumps" referred to in the previous section. The notion 
that the narrow optical lines arise in dense shocked CSM 
clumps, which is the same population that is also respon- 
sible for the lower-temperature (clumped) X-ray lines, is 
quite appealing. However the Ha lines are formed at the 
systemic velocity, while the X-ray lines are generally shifted 
by 3000-5000 km s~^ from the systemic velocity. Therefore, 
it appears that two separate populations may be needed. 



This will be investigated further in the companion paper by 
Bauer et al. (2010). 



6 DISCUSSION 

6.1 Implications for the SN progenitor 

What does our model suggest for the hydrodynamic evo- 
lution of the progenitor star? The derived density profile, 
with a low-density interior surrounded by a dense shell, and 
a wind r egion exterior t o it, resembles a wind-blown bub- 
ble ( Weaver et al.l 119771 ) formed by the interaction of two 
winds. The formation of wind-bubbles, and the evolution of 
a SN within suc h wind-bubbles, has been extensively stud- 
ied in the past (ICiotti fc D'Ercol^ll989l: IChevalier fc Lian3 
1 19891 : iTenorio-Tagle et all Il990l . Il99ll : [Franco et al.l Il991 
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iDwarkadaa |2005| . l2007al lal. These calculations have sug- 
gested that the interaction of the SN shock wave with the 
dense shell delineating the boundary of the wind bubble 
leads to increasing X-ray and radio emission. This is evi- 
dently the scenario that best resembles the emission from SN 
1996cr. This is only the second case after SN 1987A where 
a sustained increase in the X-ray and radio light curves over 
a multi-year period is seen. However, the shell in the case 
of SN 1996cr is 6 times closer to the star than the dense 
equatorial ring surrounding SN 1987A. Its density in our 
model is a few times higher than the ring density deduced 
for SN 1987A. The most recent data suggest that the shock 
has exited the shell and the emission is turning over. Thus 
SN 1996cr provides an early glimpse into what the further 
evolution of the lightcurves in SN 1987A will look like, as 
the SN shock in that case is just beginni ng to interact wit h 
the dense equatorial ring of material (R acusin et al.l 120091 ). 
In this picture, SN 1996cr represents a highly compressed, 
and much brighter version of SN 1987A. 

If we assume that the dense shell consists only of swept- 
up wind material, then the mass of the shell is equal to the 
mass of swept-up wind material. We can apply the wind 
bubble scenario to this model. In this scenario a fast super- 
sonic wind sweeps up the wind from an earlier phase into a 
thin, dense shell. We denote the earlier phase by subscript 1 
and the later phase by subscript 2. The mass of the swept- 
up shell is about 0.64 Mq. If this is made up of swept-up 
wind material, then M^sheii ~ {Mi/vi)Rsheii- Given that the 
shell outer radius is 1.5 xlO^^ cm, this gives for the exter- 
nal wind a value of Bi — Mi/vi = 8.5el5 g cm~^. If we 
assume that the wind velocity vi is 10 km s~^, as is appro- 
priate for a red supergiant wind, then the mass-loss rate of 
the star in this phase would be Mi ~ 1.33 x lO~*M0yr~^, 
which is again reasonable for a RSG wind. Thus one possi- 
bility is that the star in this phase may have been a RSG. 
If we assume a larger wind velocity, then the mass-loss 
rate must scale accordingly. Mass-loss rates on the order of 
1O~^M0 yr~^ have been deduced f or hypergiant stars with a 
wind velocity of about 40 km s~^ (|Humphrevs et al.lll997^ ■ 
and it is possible that the star in this phase could have 
been a hypergiant, consistent with the possibility that the 
progenitor was a very massive star. Even higher velocities 
would make it unreasonably large for any phase except per- 
haps t he outburst stage of a Lum inous Blue Variable (LBV) 
phase IJHumphrevs fc Davidsonlll994 V But such large mass- 
loss rates are unsustainable for a large period of time, and 
we would expect such a phase to be necessarily short-lived 
(|Humphrevs fc Davidsonlll994l l . S Doradus ty pe instabilities 
gener ally have a timescale of years to decades (|van GenderenI 
[200l|). Therefore the existence of an LBV phase would be ap- 
parent if the density profile were to change rapidly as the 
shock expanded outwards, which would show up as a strong 
deviation in the X-ray light curve from that suggested by 
us in this paper. We look forward to such measurements in 
future. 

What about the final phase of the wind before the star 
exploded? Only upper limits to the X-ray emission are avail- 
able for this period, thus only providing an upper limit to 
the density profile. In the simulations, we assume a wind 
with a value of B2 — M/vw = 6.76 xlO^^ g cm~^, which 
we refer to as the standard value of B2 . This means that for 
a wind velocity of a 1000 km s~^, the mass-loss rate is on 



the order of 10~®MQyr~^ (Figure [7)|. The wind density is 
most constrained by the upper limits on the X-ray emission 
between days 600-900. A factor of 18 increase in the value 
of B2 delays the high-flux turn-on and increases the discrep- 
ancy for the year 2000 measurements by 20 to 30 % . Thus 
the maximum mass-loss rate that would be able to produce 
a reasonable light curve (within 30% of the observed values) 
and just about satisfy other constraints (for a velocity of 
1000 km s"^) would be about 1.8 xlO~^M0yr~\ 

At the lower end we find that the density could be sig- 
nificantly smaller before the shock collides with the dense 
shell too early and exceeds the upper limits on days 700- 
900 (although this can be somewhat mitigated by changing 
the explosion date, which is uncertain by about a year). For 
the same 1000 km s^^ velocity, we find a lower-limit to the 
mass-loss rate of about 3.26 xlO~^M0yr~^. Thus, although 
we cannot exactly constrain the inner density, we find that 
the best fit light curve leads to mass-loss rates between 3. 
xlO"^M0yr"^ and 2. xlO"'^M0yr"\ for wind velocities of 
1000 kms-\ 

At the higher end, this mass-loss rate would be com- 
patible with a Wolf-Rayet star progenitor for the SN. This 
would make the SN a Type Ib/c. At the lower end, and 
especially for lower velocities, the mass-loss rate and wind 
velocity combination approaches that deduced for SN 1987A 
(Chevalier fc Dwarkada"3ll995l : lLundavist|[l999l : iDwarkadasI 
2007c b). and may further strengthen the analogy between 
SN 1996cr and SN 1987A. While we are not suggesting that 
the progenitor of SN 1996cr was a blue supergiant, it lies 
within the realm of possibility. 

One implication of these low mass-loss rates is that it is 
unlikely that the progenitor star was an LBV star. LBVs are 
highly luminous, unsta ble stars, which can undergo ep isodes 
of dramatic mass- loss (|Humphrevs fc DavidsorJll994l '). lead- 
ing to the formation of circumstellar nebulae around th e 
star. LBV nebulae are generally of order Ipc (|Weisll200ll ). 
much larger than the circumstellar nebula in our simula- 
tions. LBV stars h ave often been sugge s ted as progenitors 
for Type Iln SNe J Gal-Yam et al."2007'; 'Smith' 20081; IVinkI 
l2008l : iGal-Yam fcLeon ard 2009: Miller ctal. 2009), based 
mainly on the fact that they undergo mass-loss episodes 
with a high mass-loss rate > lO~'*M0yr~'^. LBV winds are 
also slower t han W-R winds, wit h velocities on the order of 
200 km s"^ (|Vink fc Kotakll2007l '). Assuming such a velocity 
would lower the mass-loss rates derived above by a further 
factor of 5. This would lead to mass-loss rates lying between 
about 6 xlO"^°M0yr"^ and 4 xlO~''M0yr"^ (see Figure 
7) for a presumed LBV progenitor velocity. These are sig- 
nificantly lower than mass-loss rates proposed for LBV pro- 
genitor stars. Mass-loss from an LBV in the quiescent stage 
may be consistent with the higher e nd of the rates deduced 
here (|Hum phrevs fc Davidsonlll994r ). but then it is impossi- 
ble to distin guish this from other massive supergiant stars 
(jStothers fc ChirJI 19961 ) based purely on hydrodynamic con- 
siderations. In any case by assumption this would indicate 
that the wind bubble was not the result of an LBV outburst. 

We can further constrain the wind parameters by in- 
cluding the shell formation time, t2. As above, we assume 
that the shell was formed by the interaction of two winds 
with constant mass- loss parameters, the first wind being 
the earlier RSG wind, followed by the second wind phase. 
The shell radius is given by Rsheii = {L2/[1.5Bi]Y'^t2 cm 
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(|Dwarkadaa 1 19971 ) where L2 = 0.5M2i'tu2^ is the mechani- 
cal luminosity of the second wind. With Raheii ~ 1-5 xlO^^ 
cm, and B\ = 8.5 x 10^* g cni~^, we still have one more 
unknown, the age of the shell. For an age i2 of 10,000 years, 
we obtain L2 ~ 1.37 xlO'^'^ ergs s~^. This is on the weak 
side for a W-R star, but may be comparable to that for the 
BSG progenitor of SN 1987A. Using the standard value of 
B2 above, this gives a wind velocity Vui2 of 160 km s~^, and 
mass-loss rate of the second wind about 1.7 xlO~^A/Qyr~^. 
For the lowest value of B2 we get a velocity of 1098 km s~^ 
and mass-loss rate of 3.6 xlO'^MQyr"'^, while for the high- 
est value we get a velocity of 61 km s~^ and a mass-loss 
rate of 1.17 xlO'^^Moyr'^ Values of B2 shghtly lower than 
the standard value are comparable to the BSG progenitor 
of SN 1987A, while other combinations do not seem rea- 
sonable for known stars. Thus an age of about 10^ years is 
possible only for low values of B2, leading to a low velocity 
and low mass-loss rate progenitor. The expansion velocity 
of the shell Rsheii/t2 is about 5 km s~^ (which technically 
requires that the RSG wind velocity be also reduced, but 
factors of a few are easily accommodated). 

It can be seen that L2 oc t~^, and therefore if we re- 
duce the age t2 to 1000 years, we get L2 ~ 1.37 xlO^® ergs 
s~^. This implies a mass-loss rate M2 of 1.7 xlO~®M0yr'^ 
and a wind velocity Vm2 of 1594 km s~^ for the standard 
value. The expansion velocity of the shell now is around 48 
km s~^ . These parameters are reasonable for a W-R star. It 
can be seen that lowering the values of B2 by more than a 
factor of 30 increases the wind velocity to > 5000 km s~^, 
which is incompatible with known wind velocities. The high- 
est value of B2 gives a velocity of 608 km s~^ for a mass-loss 
rate of 1.16 x lO~^M0yr~^. In this case it seems that values 
ranging from about 3% of the standard value to 4 times the 
standard value give mass-loss parameters compatible with 
known stars. 

If we reduce the age t still further, thereby increas- 
ing 1/2, the wind velocity increases proportionally. Reduc- 
tion in age by a factor of 2 gives a mass-loss rate of 3.4 
xlO~''M0yr~^ for a wind velocity of 3188 km s~^ for the 
standard value of B2. Values much lower than the standard 
value are pretty much ruled out because they give too high 
wind velocities. At the higher end this gives a wind velocity 
of 1216 km s~^ with a mass-loss rate of 2.33 x lO~^M0yr~^, 
probably reasonable for a W-R star. It is unlikely therefore 
that the age can be reduced by even a factor of 2 for the low- 
est values of B2. For high values of B2 the age could possibly 
be reduced by a factor of 3-5. An age of 250 years with the 
highest possible value of B2 gives a wind velocity of 2431 km 
s~^ with a mass-loss rate of 4.67 lO~^M0yr~^. These may 
just about be compatible with W-R stars. However we note 
that ages less than 1000 years gives velocities for the swept- 
up shell that are much larger than observed shell velocities 
around massive sta rs. Observed shell velocities of nebulae 
around W- R stars |Cappa et all 120031) . as well as around 
BSG stars (jCrotts fc Heathcotdll99ll ). are all in the range 
of 10 km s~^ or less. Alt hough muc h larger shell velocities 
are seen in LBV nebulae l|Weisll200ll ). the mass- loss rate and 
velocity combinations that are possible here never resemble 
those seen in LBV outbursts. The conclusion is that the age 
of the shell lies between 10^ and 10* years, with a preference 
for a smaller age. There is a large range over which a W-R 
progenitor seems viable, and a smaller range where a SN 



1987A-like BSG progenitor may work. The range of values 
for the mass-loss and wind velocity, given the parameters 
used and the acceptable age of the wind-blown shell, are 
shown in Figure [T] 

The above analysis then suggests the following scenario 
for the progenitor: It was a massive O star that evolved off 
the main sequence (MS) to become a RSG or hypergiant 
with a mass loss rate of ~ lO~*M0yr~^. This star then 
entered a second, brief post-MS phase wherein it became 
either a BSG star, or more likely a W-R star. This second 
post-MS phase lasted between 10^ and 10* years, with a 
preference towards the smaller age for a W-R progenitor, 
and a larger age for the 87A-like BSG progenitor. The BSG 
case suggests a star with a zero ag e main-sequence (ZAMS) 
of 15-30 Mq ijWooslev et al.lll98a ). The W-R case points to 
a massive ZAMS star > 3OM0, which went from a RSG to a 
brief W-R phase before exploding. The abundances derived 
from the X-ray fitting ( ^6.2p suggest a preference for the 
W-R case. 

A final phase lasting for upto 10* years is not unlike 
the case of SN 1987A, where the final BSG phase is as- 
sumed to hav e lasted for about 20,000 years before the 
star exploded (|McCravll2do3 . 120071 ') . The small timescale de- 
rived for the final pre-SN phase puts SN 1996cr in the cat- 
egory o f several other recent SNe, including SN 1987A, SN 
2006JC (iFolev et al."2007'; Tominaga et al. 2008), SN 2006gy 
(Smith fe McCrav 2007), SN 200 8iy (Miller et al. 2009 ) and 
SN 1994W (|Chugai et allliooi ). which lost about a solar 
mass of material a short while before the stellar explosion. 
In either case the current observations and our simulations 
are helping to shed light on the last < 10* years , or < 1% 
of the star's multi-million year history. 

Is it possible that the wind bubble scenario does not 
apply at all, that the shell is expanding ballistically, and 
that it was ejected in an LBV er uption? Giv en the speeds 
of the winds seen in LBV nebulae (|Weisll200ll ). and the fact 
that the shell is expanding ballistically, one would expect it 
to have speeds of at least 500 km s~^. The above analysis 
shows that such an eruption would have occurred less than 
a 100 years before the SN explosion. However, if the shell is 
purely wind driven, it is not clear why the interior is such low 
density, and why the shell is deduced to be so thick. Density 
contours from radiatively driven win d models of LBV nebu- 
lae can be seen in Figures 3 and 4 of iDwarkadas &: Owockil 
(2002), and do not resemble the density structure derived 
herein, although it is possible that other mass-loss mecha- 
nisms may result in different density distribution. In partic- 
ular the thickness of the shell (50% of the radius) implies 
some kind of wind interaction, and would be hard to ex- 
plain in a model where the wind was expanding ballistically. 
Furthermore, the fact that the mass in the shell is found to 
be equal to the mass in the external wind of the assumed 
density upto the outer radius of the shell, would be an ex- 
treme coincidence if the shell was not composed of swept-up 
external wind material. Therefore we consider this to be a 
less likely possibility. 



6.2 Abundances 

Table 1 describes the abundances that we have derived in 
order to compute spectra that match the observations. The 
statistical 1-cr range of values for each abundance is also 
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Figure 7. The range of mass-loss parameters for the pre-SN wind, assuming the wind density parameters Bi (denoted by an orange 
star symbol) and B2 (see i|6.ip . and timescales between 1000 and 10* years. Possible values lie at the intersection of the B2 and 42 lines. 
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green triangles. Various portions of the diagram that correspond to mass-loss parameters for W-R stars, SN 1987A-type BSGs, LBVs, 
RSG and Hypergiant stars are marked. The most likely progenitor is a W-R star. 



given. Figure [8] graphically displays the range of abundances 
deduced from our X-ray fitting. The H, He, N, O values are 
shown with very small error bars to indicate that they are 
fixed (H, He) or undetermined (N,0). In the figure we also 
compare the abundances to those derived from LETG and 
HETG observations of SN 1987A JZhekov et al.ll2009l'). as 
well a s those given for WCE stars bv lMellema &: Lundgvistl 
1(2003) ■ These comparisons lead to the following observa- 
tions: (1) The X-ray data provide very little constraint on 
the ratios of He to H, or /i, of the plasmas. Agreement with 
the light curve and spectra was improved with our choice of 
a H-rich GSM and a more He-rich ejecta. However it is likely 
that other possibilities, e.g., both plasmas with a /x ~ 1.0, 
could be accommodated with some adjustments to the den- 
sity profile. (2) Our values for the relative abundances of the 
heavier elements in the GSM are consistent with those for 
the GSM around S N 1987A derived fro m fitting the LETG 
and HETG spectra l|Zhekov et al.ll2009l '). Those abundances 
were found to be consistent with those of the LMG, and 
with SNRs in the LMG. In the case of the Gircinus galaxy, 
lOliva et al.l (|l99 9) arrive at an average metallicity of -0.7 ± 
0.3 from considerations of the equivalent width of the stel- 
lar GO lines in a square of lOOpc on a side. This indicates 
a metallicity lower than that of the LMG, which would en- 
hance the potential for a BSG to explode. (3) The main 
difference between SN 1996cr and SN 1987A abundances is 
in the H and He that make up the "continuum". Our fits 
predict H and He abundances that are an order of magni- 



tude lower than those derived for SN 1987A. (4) The low H 
and He in the ejecta are suggestive of a W-R star which has 
expelled its H and He envelope, consistent with the analy- 
sis from purely hydrodynamical considerations above. The 
abundance measurements are close s t to th ose for a WGe star 
as given in iMellema fc Lundgvistl l|2002l 'l (5) The GSM in 
this case, although mainly composed of H and He, is over- 
abundant in heavy elements like S and Si. The existence 
of strong Si and S H-like lines, and the very strong Fe H- 
like line at ~ 7 KeV in the 2001 XMM spectrum, requires 
a metal-rich plasma. Although RSGs are not expected to 
have such metal-rich winds, it is interesting to note that the 
RSG WOH G64 shows the presen ce of [SII] 6717/6731 in 
its spectrum IjLevesque et al.ll2009l ). while the M supergiant 
VY GMa shows the presenc e of Si I li nes, and SiO maser 
emission (Wallcrstcin & Go nzalej|200ll ). along with various 
signatures of other metal species. Therefore one possibility 
is that the winds are enhanced in Si and S. An alternate pos- 
sibility is that the shocked GSM could actually be a mixture 
of GSM and ejecta material, due to mixing at the interface 
between the shocked GSM and shocked ejecta, which can 
lead to Rayleigh- Taylor "fingers" of ejecta propagating into 
the GSM. In order to test this possibility, we have assumed 
a modified GSM composition, that i s made up of equal parts 
(by number) of the SN 1987A GSM jZhekov et al .' 2009) and 
SN 1996cr ejecta. Such a 50-50 mixture has a resulting abun- 
dance similar to that of our "GSM", and, if used in place of 
our GSM abundance values, gives a similar light curve and 
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spectra. Thus what we are labehng as "CSM" could be con- 
taminated by ejecta mixing, and the derived abundance re- 
flects this mixture. (6) The H abundance is somewhat higher 
in the CSM while the He abundance is higher in the ejecta. 
This may indicate that the star gave off most of its H layer 
but retained somewhat more of its He layer, pointing per- 
haps to a Type lb SN. The characteristic velo city computed 
in j}3] is also closest to that of a Type lb SN (iMaurer et al.l 
120091 ). The overall abundance of metals is still high in the 
CSM, suggesting perhaps that most of the H and He was 
lost in a much earlier stage, and is much further out in ra- 
dius. This is consistent with the fact that the circumstellar 
bubble and dense shell with which the SN shock wave was 
interacting in its first decade was emitted in the very last 
stages of the stellar lifetime, and that much more CSM ma- 
terial that was released earlier is to be found further out at 
a larger radius. If so the SN should probe this material over 
the next few decades. (7) One aspect that we have not con- 
sidered, and which may play a role in this scenario, is grain 
chemistry. RSG stars are known to have dusty winds, and 
stars with higher mass-loss rates are more likely to show the 
presence of dust (Jones 200l|). It is possible that the dense 
shell, if composed of RSG wind material, has a significant 
dust component. Under equilibrium conditions, metals can 
condense in dust grains, thus depleting certain species and 
affecting the X-ray plasma abundances. Alternatively, given 
the presence of the high density shell and the hot gas, which 
can provide ions with an energy of several keV, it is also 
possible that sputtering of the dust by hot gas as the shock 
travels through the dense shell is important. This may alter 
the X-ray plasma composition. 



T SUMMARY 

We have computed a hydrodynamic model that is able to 
explain the X-ray light curves and spectra of the young su- 
pernova 1996cr. Our intent was not to completely explore 
the parameter space - indeed, this is practically impossible 
- but to evaluate the best possible model consistent with 
all of the available data. Our model incorporates a dense 
shell of material that lies about O.OSpc from the star and 
extends radially for about 0.015pc, with a density of about 
8 X 10* amu/cc. The density interior to the shell is quite 
low, although it is not well constrained by the observations. 
The density external to this is well approximated by a wind 
profile at present. The interaction of the SN ejecta with this 
shell and wind is able to reproduce the increasing X-ray light 
curves that have distinguished this supernova, and are also 
consistent with the behavior of the radio light curves. 

We have computed the X-ray emission from our simula- 
tions, using non-equilibrium ionization conditions, and tak- 
ing into account the inequality between the electron and ion 
temperatures. Our method is both quick and accurat4j and 
allows us to compute simulated spectra at every epoch, that 
can be directly compared with the observations. We have 
compared both the full spectrum as well as specific lines. 
This is better described in the companion paper (Bauer et 

* When referring to the accuracy of the flux obtained from spher- 
ically symmetric simulations of a SN which is still unresolved at 
all wavelengths, we mean to better than about 30%. 



al. 2010, in preparation). The fiux at various epochs, and 
the computed X-ray spectra, agree well with the observed 
data, thereby validating our model. 

Our hydrodynamic model suggests that the progenitor 
of SN 1996cr was more likely a BSG or W-R star, with a fast, 
low density wind. The interaction of this wind with a prior 
wind, possibly from a RSG phase, gave rise to the wind- 
blown shell with which the SN shock wave was colliding 
in its first decade, leading to the increasing radio and X- 
ray enrission. The sinrilarity to SN 1987A in nrany different 
aspects is quite striking. An LBV progenitor is inconsistent 
with much of the data. 

Given the similarity with the medium around SN 
1987A, which has a clear ring structure surrounding it, it 
is interesting to ask whether a ring structure, rather than a 
fully developed spherical shell, would be a possible scenario 
for SN 1996cr. In our study, we have found that a somewhat 
denser shell (about 50% higher density), coupled with the 
assumption that the emission was arising from a region of 
27r steradians, would give approximately the same level of 
fiux. Thus it is possible that a donut-shaped structure rather 
than a spherical shell could be accommodated, although we 
have not investigated the detailed spectra. It seems unlikely 
that a denser shell would fit because then the VLBI radius 
constraint would not be satisfied. Further investigations in 
this direction are left to a future paper as they would require 
far more data for comparison than is currently available. 

It is possible that the further expansion of the shock 
may not agree with our supposition that the nredium outside 
the shell is a RSG wind, or indeed any kind of wind. Only 
further data on the expansion can shed light on this, and 
we eagerly look forward to continually exploring the X-ray 
emission of this interesting and enigmatic object with the 
Chandra and XMM-Newton telescopes e.g. through yearly 
fiux monitoring and and future HETG spectroscopy if feasi- 
ble. Similarly, a second VLBI measurement, which will pro- 
vide a second radius measurement and a direct indication 
of the current shock velocity, is eagerly awaited. Finally, we 
also have applied for further optical and IR observations 
that will help to create a nrulti-wavelength history of this 
SN, and enable us to understand the various properties of 
this fascinating object. 

The relationship of SNe to their progenitor stars is quite 
uncertain, as mentioned in the introduction. Surveys that 
try to pin down the progenitors are doing a fantastic job 
fsee ISmarttll2009l . and references within), but the problem 
is difficult. Theoretical considerations like those used in this 
paper provide another completely different technique to ex- 
plore the CSM around the supernova, and link it to its pro- 
genitor star. No doubt there is uncertainty in this, as in any 
other method that is trying to ascertain the properties of 
a star that no longer exists. But by using the same hydro- 
dynamic simulation to fit multi-epoch X-ray spectra as well 
as light curves without adding too many free parameters, 
we have tried to mitigate uncertainties as far as we can. 
The fact that the simulations reproduce the observed shape 
of the continuum at each epoch without any additional pa- 
rameter speaks for the validity of our model. It is only by 
amalgamating the results from various methods that we can 
hope to understand the true progenitors of the various types 
of SNe. In the case of a SN such as 1996cr, where the possi- 
bility of looking for a progenitor star at the same location. 
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Table 1. The ejecta and CSM abundance distributions as deduced from spec- 
tral fitting. Th e abundances are ex p ressed relative to Si with respect to the 
values listed in I Anders fc Grevessa 1 119891) : the 1-a range for elements with 
detectable lines are also given. The "continuum" abundances, H through O, 
are not well constrained by the data (see text). Other columns give these 
abundances converted to number fraction and to percent by mass. For the 
ejecta (CSM) the mean ion mass, ^ia, is 3.548 (1.693) amu and the mean 
mass per particle, /i, is 1.236 (0.755) amu. 









Abund.s w.r 


t. Si (AG89) 


Number fraction 


Percent mass 






Element 


Ejecta 


CSM 


Ejecta 


CSM 


Ejecta CSM 




H 


0.0121 


0.0811 


2.02C-01 


7.96e-01 


5.69 47.0 






He 


0.483 


0.207 


7.89e-01 


1.98e-01 


90.0 46.8 






C 


0.128 


0.336 


7.77C-04 


1.20e-03 


0.263 0.849 






N 


0.128 


0.336 


2.40C-04 


3.70e-04 


0.095 0.306 






O 


0.128 


0.336 


1.82C-03 


2.81e-03 


0.821 2.655 




Ne 


2.09 [1.74-2.48] 


0.07 [0-0.68] 


4.30e-03 


8.81e-05 


2.427 0.104 






Mg 


0.82 [0.70-0.96] 


0.46 [0.13-0.85] 


5.23e-04 


1.70e-04 


0.354 0.242 






Si 


1.00 [0.91-1.10] 


1.00 [0.81-1.26] 


5.93e-04 


3.48e-04 


0.468 0.576 






S 


0.92 [0.77-1.08] 


1.43 [1.13-1.78] 


2.49e-04 


2.28e-04 


0.225 0.431 






Ar 


0.82 [0.51-1.15] 


1.22 [0-2.03] 


5.00e-05 


4.34e-05 


0.056 0.102 






Ca 


0.13 [0-0.48] 


2.82 [0-3.94] 


5.04e-06 


6.35e-05 


0.006 0.150 






Fe 


0.49 [0.44-0.55] 


0.46 [0.41-0.52] 


3.83e-04 


2.10e-04 


0.605 0.693 
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Figure 8. Abundances derived from our X-ray fitting, compared to others found in the literature. The H, He, N, O values are shown 
with very small err o r bars to indicate that they are fixed (H, He) or undetermined (N,0). The Red hexagons are SN 1987A measurements 
from IZhekov et al.l ||2P03). SN 1996cr appears to have 10 times more Z > 8 metals compared to the He (and H) there, and seems to 
posses a H, He-poor or metal-rich CSM compared to SN 1987A. But the ratio of the heavy (Z > 8) metals is about the same - the Oxygen 
value in SN 1996cr for instance could be easily set equal to that seen in SN 1987A without affecting the light-curve fits. 



which faded away after the SN explosion, no longer exists, a 
theoretical analysis may be the only method of excavating 
the SN progenitor. 
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8 APPENDIX - CALCULATION OF THE 
PROGENITOR WIND PARAMETERS 

A detailed description of the wind parameters in an inter- 
acting winds scenario , for various shock stre n gths a nd other 
variables, is given in IChevalier fc Imamural 1 19831 ). A sim- 
plified description, as outlined below, is sufficient for our 
purposes. 

We start with the equation for the radius of the shell in 



Rsheii = {L2/[1.5Bi]y^h2 (9) 

Noting that L2 ~ 0.5M2Vw2 s^nd «„2 ~ M2/B2 we get 

M2t2 



R 



B2 



which gives 



{M2t2Y 



(10) 



(11) 



Note that M2t2 is the total mass of the emitted wind 
over time t2 . This includes the mass of the freely expanding 
wind Mfyj, which extends upto the wind termination shock 
Rt {Mfw — B2*Rt) and the mass of the shocked wind, which 
stretches from Rt to the inner boundary of the shell Ri. If 
we assume that the density of the shocked wind is constant 
throughout, a reasonably good approximation, and that the 
wind termination shock is strong, i.e. with a shock jump of 
4, then the density of the shocked wind is 4 times that of 
the free wind at Rt, psw ~ B2/{4,-kR^) * 4, and the mass of 
shocked wind 

B2 47r, 



M, 



i^R't) 



^—{Rl-Rt 



(12) 



Therefore, with M2t2 
B2 



Mfw + Msw we get an equation for 



'-i 



R^ 



B2[Rtil+l{^. 



m' 



sRLBi 







(13) 



This equation can be solved using standard techniques 
for all Rt between and Ri to get appropriate values of B2. 

Our calculations constrain the density of the shocked 
wind Paw, so we compute the values of Rt and B2 that pro- 
vide us with the constrained values of psw, and thus set 
limits on the progenitor wind density. 
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